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ABSTRACT 



This is the first paper in a series dedicated to the study of the emission-line and 
CN I continuum properties of narrow hne Seyfert 1 galaxies (NLSls). We carried out a 

systematic search for NLSls from objects assigned as "QSOs" or "galaxies" in the 
K>" I spectroscopic sample of the Sloan Digital Sky Survey Data Release 3 (SDSS DR3) 

ly-^ I by a careful modeling of their emission lines and continua. The result is a uniform 

I sample comprising ~ 2 000 NLSl. This sample dramatically increases the number 

O ' of known NLSl by a factor of ~ 10 over previous compilations. This paper 

O 



presents the parameters of the prominent emission lines and continua, which 
were measured accurately with typical uncertainties < 10%. Taking advantage 



Qh| of such an unprecedented large and uniform sample with accurately measured 

O ' spectral parameters, we carried out various statistical analysis, some of which 

' 

^ ' were only possible for the first time. The main results found follows. 

^ ' (1) Within the overall Seyfert 1 population, the incidence of NLSls is strongly 

dependent on the optical. X-ray, and radio luminosities as well as the radio- 
^ ' loudness. The fraction of NLSl peaks around SDSS g-band absolute magnitude 

' Mg ~ —22™ in the optical and ~ 10"^^'^ erg in the soft X-ray band, and 

decreases quickly as the radio-loudness increases. (2) On average the relative 
Fell emission, 7^4570 = Fe//(4434 - 4684) /H/3, in NLSls is about twice that 
in normal AGNs, and is anti-correlated with the broad component width of the 
Balmer emission lines. (3) The well-known anti-correlation between the width of 
broad low-ionization lines and the soft X-ray spectral slope for broad line AGNs 
extends down to FWHM ~ 1 000 km in NLSls, but the trend appears to 
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reverse at still smaller line widths. (4) The equivalent width of H/3 and Fell 
emission lines are strongly correlated with the H/9 and continuum luminosities. 
(5) We do not find any difference between NLSls and normal AGNs in regard to 
the Narrow Line Region. (6) We have examined the black hole mass vs. stellar 
velocity dispersion {Mbh — <7*) relation for a subsample of 308 NLSls for which 
cr* could be measured directly from fitting the starlight in the SDSS spectra with 
our stellar spectral templates. A significant correlation between Mbh and cr* is 
found, but with the bulk of black hole masses faUing below the values expected 
from the Mbh — cr* relation for normal galaxies and normal AGNs. This result 
indicates that NLSls are underage AGNs, where the growth of the SMBH lags 
behind the formation of the galactic bulge. (7) We also find that the FWHM 
of [Nil] line is well correlated with a* in 206 NLSls, for which both parameters 
could be derived with reasonable accuracy. The [Nil] width can predict the stellar 
velocity dispersion to an accuracy of ~ 30%. A similar Mbh — o"* relation could 
be found for a larger sample of 613 NLSls on making use of the more reliable 
measurements of FWHM[NII]. 

Subject headings: galaxies: active — galaxies: Seyfert - galaxies: emission lines 

1. Introduction 

Active galactic nuclei (AGNs) — Seyfert galaxies and quasars — are characterized by non- 
thermal continuum emission over almost the whole electromagnetic spectrum from their 
nuclei and prominent emission lines in the ultraviolet and the optical. The overall spectral 
energy distribution (SED) can be described as an underlying power law of the form oc 
interspersed with bumps and dips. The AGNs can be classified into two types according to 
their emission spectrum: type 1 AGNs show both broad emission lines with full width at half 
maximum (FWHM) of a few thousand km and narrow emission lines of a few hundred 
km s~^, while type 2 AGNs show only narrow emission lines. Permitted and semi-forbidden 
lines are seen with both broad and narrow profiles, while forbidden lines are seen only with 
narrow profiles. It was recognized at the end of 1970s (e.g., Osterbrock 1978) that such a 
split might be an orientation effect involving obscuration of the broad- line region (BLR). 
According to unification models, all AGNs have intrinsically the same physical structure, 
but in type 2s the BLR and continuum emission region are obscured by a presumed dusty 
torus located somewhere between the BLR and the narrow line region (NLR). 

It seems that the division of AGNs into type 1 and type 2 is well-defined and the 
evidence for unification schemes is compeUing (cf, Antonucci 1993). Interestingly, almost at 
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the same time when the unified models were first brought forward, it was found that two 
objects, Mrk 359 (Davidson & Kinman 1978) and Mrk 42 (Koski 1978; PhiUips 1978), have 
all the spectral properties of a typical type 1 AGN, except that the "broad permitted lines" 
are comparable in width to the forbidden lines as in the case of type 2 AGNs. Some other 
unusually narrow "broad line emitters" (e.g., Mrk 493 and Mrk 783, Osterbrock & Dahari 
1983) were identified in the following years. This intriguing mixture of features makes this 
"pecuhar" group of AGNs very interesting. Osterbrock & Pogge (1985, hereafter OP85) 
denoted these objects as narrow fine Seyfert 1 galaxies (NLSls). Their original classification 
included the following two empirical and somewhat subjective criteria: 1) the Balmer lines 
are only shghtly broader than forbidden fines; and 2) the line ratio [OIII]X5007/HP < 3; 
but exceptions are allowed when strong high- ionization iron lines, such as [FeVII]A6087 and 
[FeX]A6375 are present, which are seldom seen in type 2 AGNs. In his spectropolarimetric 
study of NLSls, Goodrich (1989) treated them as an extension of normal broad-line Seyfert Is 
(hereafter BLSls) to lower FWHM of permitted lines and quantified the first criterion of 
OP85 with FWHM{H[3) < 2000 km s'^. 

Apart from the unusually narrow Balmer lines, NLSls often show the strong permitted 
Fell emission lines in their optical and ultraviolet spectra; they also show some other extreme 
properties, such as a steep soft X-ray spectrum, rapid X-ray variability, dense line-emitting 
gas, and a commonly blue shifted UV line profile (Boiler et al. 1996; Leighly 1999; Aoki & 
Yoshida 1999; Wills et al. 1999; Leighly 2000). In fact, these properties are extensions of 
the following correlations between various observables and the if/3 line width (Boroson & 
Green 1992; Laor et al. 1994; Wang, Brinkmann, & Bergeron 1996; Xu et al. 2003): the 
narrower the Hfi line, the steeper the X-ray spectrum, the stronger the Hf3 blue wings, 
the faster the X-ray variation, the stronger the optical Fell emission, the weaker the [0111] 
lines, and the denser the BLR clouds. These form the so-called eigenvector 1 (El, Boroson 
& Green 1992; Sulcntic ct al. 2000). El is taken as the strongest and most far-reaching 
type 1 AGN unification yet discovered (e.g., Marziani et al. 2001), which incorporates the 
geometry, kinematics, and physical condition of the central regions in AGNs. Its underlying 
physical driver is related to either the basic physical parameters of the black hole accretion 
disk such as the black hole mass and accretion rate, or geometrical parameters such as the 
viewing angle (Boroson & Green 1992; Boiler et al. 1996; Wang, Brinkmann, & Bergeron 
1996; Laor et al. 1997; Sulentic et al. 2000; Marziani et al. 2001). Since NLSls lie at the 
extreme negative end of El, they may well occupy some extreme regions in the parameter 
space of some primary physical parameters, whatever they are. We can, therefore, use the 
unusualness of NLSls to test the viability of the different AGN models. 

However, some of the basic properties of NLSls remain unclear even though it has been 
two decades since their first identification. We summarize these issues under the following 
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two headings. 

Properties extensively explored, but no consensus reached yet: 1) For Seyfert 1 galax- 
ies and QSOs, a correlation between the soft X-ray photon index (from ROSAT ob- 
servations) and the Hf3 line FWHM has now been well established (Laor et al. 1994; 
Boiler et al. 1996; Wang, Brinkmann, & Bergeron 1996). In comparison, for the NLSls, 
which have on average steeper photon indices, the same correlation shows a larger dis- 
persion. It is still not known whether the flatness of the X-ray spectra of some NLSls 
is intrinsic or is due to internal absorption. 2) Boroson & Green (1992), Mcintosh et al. 
(1999) and Grupe et al. (1999) found that the strength of the Fell and [OIII] hnes are 
anti-correlated, but this was not confirmed by Veron-Cetty, Veron, & Goncalves (2001, 
hereafter VVGOl). It is unclear whether this correlation should be included in El. 3) 
Rodriguez-Ardila et al. (2000) claimed that the emission line ratio of the narrow line 
region (NLR), [0111] /H/9, is in the range of ~ 1 — 5 for NLSls, significantly less than 
the values ~ 10, typical of normal AGNs, while VVGOl did not find any difference 
between the two types. 

Virgin territories: Compared to the well studied emission line and X-ray properties of 
NLSls, little or no attention has been given to 1) their broad-band continuum prop- 
erties; 2) the relationship between the NLSl phenomena and the evolution of the host 
galaxies; 3) the type 2 counterparts of NLSls predicted by AGN unification models; 
and 4) their spectral variability. 

To resolve the controversies and to seek new insights, it would be helpful to have a large- 
size homogeneous sample of optically selected NLSls with accurately measured parameters 
of both narrow and broad lines. The Sloan Digital Sky Survey (SDSS, York et al. 2000) 
is expected to be able to provide such NLSl samples resulting from well-defined selection 
criteria that are large enough for serious statistical studies; furthermore, it is expected that 
potentially peculiar and interesting objects will be identified and studied with follow-up 
observations. The first attempt in this direction was made by Williams, Pogge, & Mathur 
(2002, hereafter WPM02) who compiled a sample of NLSls from the SDSS Early Data 
Release (EDR, Stoughton et al. 2002), which comprised 150 objects. With a roughly 10- 
fold increase in size over the EDR, the presently available SDSS spectroscopic data can be 
expected to yield many more NLSls for the study. 

We carried out a systematic search for NLSls from the spectroscopic sample of the Sloan 
Digital Sky Survey Data Release 3 (SDSS DR3, Abazajian et al. 2005). In order to obtain 
from SDSS spectra accurate and reliable fine widths that are essential for the selection and 
the follow-on analysis, we carefully fitted emission lines, AGN continuum and starlight in 
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a self-consistent manner. As a result, we now have a uniform sample comprising as many 
as 2 Oil NLSls. The accurately measured emission-line parameters and hence the selection 
of sample objects are an improvement on WPM02, rendering our sample less affected by 
selection biases compared to WPM02. This paper is the first of a series which reports on our 
ongoing project for extensive studies of NLSls by fully exploring the potential of this large- 
size sample. Here we present the object catalog, which will be kept updated as new SDSS data 
are further available; also presented are some preliminary statistical results on the sample as 
a whole. We shall attempt to address some of the issues raised above in forthcoming papers. 
The plan of this paper is as follows: In the next section, we describe in detail the methods 
we used to process all the low redshift QSOs and galaxies in the SDSS DR3. The sample 
selection is described in §3. The properties of the sample arc discussed in §4, including a 
statistical analysis of the emission lines and remarks on potentially interesting individuals. 
We explore the black hole-bulge relation for NLSls in §5. Our main results are summarized 
in the last section, followed by prospects of future work. A A-dominated cosmology with 
Hq = 70 km Mpc~^, Qm = 0.3 and Qa = 0.7 is assumed throughout the paper. 



2. Optical Spectral Analysis 
2.1. An overview 

Our NLSls were selected from the SDSS DR3 catalogs that are classified as "galaxy" 
or "QSO" by the SDSS spectroscopic pipeline. Only objects with a redshift z < 0.8 were 
considered to ensure that the HP and [OIII] emission fines would fall within the wavelength 
coverage of the SDSS spectrograph (Details of SDSS spectroscopy can be found in Stoughton 
et al. 2002). The spectra of the resulting 387 483 extragalactic objects (12 824 "QSOs" and 
374 659 "galaxies") were first corrected for the Galactic extinction using the extinction 
curve of Schlegel et al. (1998), and then shifted back to their rest frames using the redshifts 
provided by the SDSS pipeline. 

Contamination from stellar light: The fibers that feed the SDSS spectrograph subtend 
a diameter of 3" on the sky, which corresponds to ~ 6.5 kpc at z = 0.1. AGN spectra taken 
through such a fixed aperture may include significant light from the host galaxies, as exhib- 
ited in the composite spectra of SDSS quasars (Vanden Berk et al. 2001). Careful removal 
of the starhght is essential for reliable measurement of the emission fines, and hence correct 
identifications of NLSls from the huge SDSS data. Actually, the starhght "contamination" 
provides valuable information about the host galaxies, which is of interest in itself. For this 
and other purposes, we have developed a technique of properly modeling the stellar com- 
ponent (see Lu et al. 2006 for details of the method). This method decomposes the SDSS 
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spectra of active galaxies into stellar and non-stellar nuclear components, provided that the 
two components are comparable in strength. A description of this method is given in §2.2. 

Broad-line and narrow-line decomposition: The Balmer emission lines of NLSls may 
include significant contribution from narrow- line components (WGOl; Rodriguez- Ardila 
et al. 2000). In the selection of NLSls based on the hue's FWHM, it is therefore more 
appropriate to use the broad component rather than the total emission line, as was done in 
WPM02. Moreover, the dc-blended narrow components have interest of their own for the 
study of the NLR properties in NLSls. We therefore performed emission-line decomposition 
in the analysis. For reliable line decomposition, good spectral quality is required in terms of 
both S/N ratio and spectral resolution. Rodn'gucz-Ardila et al. (2000) demonstrated that 
spectra with a high S/N ratio (the exact value is unclear but likely to be around 30/pixel) and 
a moderate resolution (340 km FWHM at Ha) are sufficient for the decomposition of the 
permitted emission-lines in NLSls. In comparison, the SDSS spectra used in this work have 
a somewhat lower S/N on average but better resolution (estimated to be 132 ± 12 km 
FWHM at Ha, cf. Greene & Ho 2005); the effect of the latter compensates the effect of 
the former. Spectral decomposition was also carried out in VVGOl, in which the spectral 
resolution {FWHM ~ 3.4 A) was actually a little bit worse than that of the SDSS spectra. 
Wc therefore anticipate that the spectral quality of the SDSS data is adequate for our 
purposes, at least for the majority of objects. A detailed account of spectral-line fitting is 
presented in § 2.3. 

Spectral fitting procedure: In the processing of this large data set wc employed a two-step 
iteration procedure, which is outlined below and detailed in the following two subsections. 
The continua of the broad emission-line AGN candidates were first decomposed into starlight 
and nuclear components. Beforehand, two kinds of spectral regions had been masked out: 
firstly, bad pixels as fiagged by the SDSS pipehne; secondly, the wavelength ranges which 
may be seriously affected by prominent emission lines characteristic of QSO spectra. To 
determine the latter, the composite SDSS QSO spectrum (Vandcn Berk et al. 2001) was 
used during the first iteration. A pseudo-continuum was fitted and subtracted, which was 
chosen as a non-negative linear combination of properly generated templates taking into 
account the contribution from both starlight and nucleus (§2.2). The leftover emission-line 
spectrum was fitted with Gaussian and/or Lorentzian profiles, separating broad and narrow 
line components (§2.3). Then the measured emission line spectrum was used to replace the 
composite SDSS spectrum in the second iteration. This procedure was reiterated until the 
fitted parameters of both the continuum and the emission-line spectrum converged to an 
acceptable accuracy. 
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2.2. Decomposing Starlight and Nuclear Continuum 

We modeled the pseudo-continuum of the 387 483 SDSS galaxies and QSOs with 2 to 4 
components 

S{X) = Anost{E'^%, A) A{X) + Awe«.(^r-t"', A) [bB{X) + cC(A) + dD{X)] (1) 

where S{X) is the observed spectrum. A{X) = Yl^=i^i ^Ci{X,a^) represents the starlight 
component modeled by our 6 synthesized galaxy templates, which had been built up from 
the spectral template hbrary of Simple Stellar Populations (SSPs) of Bruzual & Chariot 
(2003, hereafter, BC03) using our new method based on the Ensemble Learning Independent 
Component Analysis (EL-ICA^) algorithm. The merits of these galaxy templates are that, 
firstly, their physical meanings are clear; secondly, in most cases, an over-fit to the spectra 
can be avoided because the most prominent features of the BC03 SSP library are embodied 
to the greatest extent in only 6 non- negative independent components (ICs). The details 
of the galaxy templates and their applications are presented in Lu et al. (2006). A{X) was 
broadened by convolving with a Gaussian of width a* to match the stellar velocity dispersion 
of the host galaxy. The un-reddened nuclear continuum is assumed to be B{X) = X~^^^ as 
given in Francis (1996). C(A) denotes the optical Fell templates of Veron-Cetty et al. (2004) 
with the width fixed to that of the broad component of Hp. D{X) represents the templates 
of higher-order Balmer emission lines (10 <n < 50) and Balmer continuum generated in the 
same way as in Dietrich et al. (2003). Ahost{E'^"Iv , A) and Anude«s(-E'^"^y"^, A) are the color 
excesses due to the extinction of the host galaxy and the dusty torus, respectively, assuming 
the extinction curve for the Small Magellanic Cloud (SMC) of Pei (1992). The fitting was 
performed by minimizing the with E^°fy, E^^y'^^, ai, cr*, b, c, and d being non- negative 
free parameters. 

To reduce the amount of computation, E'^°_^y, E^_^y^'^ , c and d were fixed at zero during 
the first iteration. Then the modeled continuum was subtracted from the observed spectrum 
and the emission line parameters were measured in the way described in §2.3. E^fy and 
j^micieus f]~gg parameters when the continuum was fitted the second time. 

For those objects in which the broad emission fine is rehably detected (at the > 5 cr 
confidence level) , c and d were also set free, and the width of the Fell multiplets was fixed 
to that of the broad component. This time the wavelength ranges of the emission 



^ICA is a new Blind Source Separation (BSS) method for finding underlying statistically independent, 

non-Gaussian components from multidimensional data. This statistical and computational technique can 
be regarded as an extension to Principle Component Analysis (PCA) but is much more powerful than the 
latter. The Ensemble Learning algorithm is an approach to find an analytic approximation applicable to the 
ICA model (Miskin & MacKay 2001). 
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lines to be masked out were determined using the measured line parameters. The iteration 
continued until the fitted parameters of both the modeled stellar spectrum and the emission- 
line spectrum converged to an acceptable accuracy. Representative examples of the EL-ICA 
starlight-nucleus decomposition are shown in Figure 1. 

2.3. Emission-Line Fitting 

The emission-line fitting procedure adopted here is similar to that used in Dong et al. 

(2005), but with some modifications in the codes to improve the accuracy of the measurement 
of the broad components of and Ha, as described in the following. The emission lines 
are decomposed into a broad and a narrow component wherever the latter contributes non- 
negligibly. The [Nil] doublet are separated from Ha. 

All narrow emission lines, except the [OIII] doublet, were fitted with a single Gaussian 
because these lines are often quite weak in NLSls and the S/N is not high for most objects. 
The [0111] doublet were fitted cither with a single Gaussian if the profile is symmetrical or of 
low S/N, or with a double Gaussian otherwise. The fiux ratios of the [0111] and [Nil] doublets 
were fixed at their theoretical values. The profile and redshift of the Ha narrow component 
were assumed to be the same as those of the [Nil] and [SII] doublets; — this approach is often 
adopted in separating the narrow and broad components in the Ha-|-[NII]-|-[SII] regime (e.g.. 
Ho et al. 1997; VVGOl; and Dong et al. 2005). This assumption can be tested using type 2 
AGNs, where narrow emission lines can be measured much more easily. Here wc carried 
out such a test using a sample of ~ 3 000 type 2 AGNs with high S/N spectra from SDSS. 
Figure 2 shows FWHM[NII] (upper left panel) and FWHM[SII] (upper right panel) against 
FWHM(Ha). It can be seen that both FWHM[NII] and FWHM[SII] are statistically well 
consistent with FWHM(Ha). For this type 2 AGN sample, then, we made a second fit after 
tying together the profiles and redshifts of Ha, [Nil] AA6548, 6583, and [SII]AA6716, 6731. The 
fitted fine width was compared with that of H/3, and a good agreement was found (Figure 2 
lower left panel). FWHM(H/3) was also found to be well correlated with FWHM[OIII], 
though the scatter is much larger than in the FWHM(H/5)-FWHM[NII] correlation. We 
therefore consider it a good approximation to tie the profile and redshift of the narrow 
component of H/9 — in the case of weak lines, either narrow or broad — to the fitted values of 
Ha+[NII] or [OIII]. We refer to these schemes as "Ha to H/3" and "[OIII] to H/3", respectively. 

In the first place, each of the forbidden and permitted lines was fitted with a single 
Gaussian. Objects were assigned as candidates of broad emission-line AGNs if their Balmer 
emission line(s) were significantly broader than the forbidden lines or their Fell multiplets 
were detected at the ^ 5 cr confidence level. The spectra of more than 6 000 candidates from 
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the SDSS "galaxy" catalog were visually inspected and 1 926 turned out to be broad line 
AGNs. These, together with 12 824 objects with z < 0.8 in the SDSS "QSO" catalog, were 
carried over into the second iteration stage. After the continuum was subtracted, the Ha 
and H/9 emission lines were fitted with two Gaussians, one for the narrow and one for the 
broad component. As mentioned above, the width and redshift of the narrow component 
of Ha were tied to those of the [Nil] and [SII] doublets; and for the narrow component of 
H/3, the "Ha to H/?" (for z < 0.39) or "[OIII] to H/5" (for 0.39 < z < 0.8) schemes were 
adopted. The objects were reserved whose broad component(s) of Ha or H/3 were detected 
at the > 10 (J level with the FWHM less than 3 000 km s~^. All the rejected objects were 
double-checked for confirmation by visual inspection to ensure that no true NLSls were 
missed (actually such checks were made throughout the whole selection process whenever an 
object was to be rejected). We also rejected at this stage the narrow emission line galaxies 
that were mis-classified as "QSOs" by the SDSS pipeline. At the stage of the third iteration, 
the broad components of Ha and Hf] were fitted with a Lorentzian profile, while the narrow 
lines were fitted in the same way as in the second iteration. Only those having FWHM(Ha) 
or FWHM(H/3) less than 2 500kms~^ were carried over into the last iteration stage. The 
emission lines were then fitted using one of the following four schemes in order of decreasing 
preference: 1) the broad and narrow components of H/3 are de-blended freely; 2) the "Ha to 
H/3" scheme; 3) the "[OIII] to H/?" scheme; and 4) Ha and H/? lines are fitted with a single 
Lorentzian (for objects having no isolated narrow emission line detected or having too noisy 
spectra). The fitting results were visually inspected before acceptance. Only 42 objects 
could not be fitted satisfactorily by using the above procedure. Among those exceptions, 
the majority have the broad components of Ha or H/3 too narrow to separate from the 
narrow components; the remaining have their redshifts incorrectly assigned by the SDSS 
pipeline. Their spectra were fitted manually. Figure 3 shows a few examples of emission-line 
decomposition in the Ha and H/3 wavelength ranges. 

3. Sample Selection 

3.1. Selection Criteria 

Based on the emission line parameters measured, we compiled our NLSl sample. The 
only criterion we adopted for identifying a NLSl galaxy is: the "broad" component of H/? 
or Ha is detected (at the > 10 a confidence level) and is narrower than 2 200 km in 
FWHM. This criterion resulted in 2 Oil NLSl candidates, whose main spectral parameters 
are summarized in Table 1. 

It is noticeable that our selection criterion is somewhat different from the "conven- 
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tional" NLSl definition (see the Introduction) in two aspects: firstly, the hne width cut-off 
of 2 200 km is slightly higher than the commonly used 2 000 km s~^; secondly, we did 
not restrict ourselves to objects with [OIII]/ H j3 < 3. These modifications were made based 
on the following considerations: (1) The line width distribution for the Balmer lines reveals 
a rather smooth transition between NLSls and BLSls, which means that any dividing line 
is somewhat arbitrary and 'operational'. Besides, different line profiles were used by the 
different authors (see WGOl and the references therein), and sometimes the narrow com- 
ponent was not taken into account properly (e.g., WPM02), resulting in underestimation of 
the hne width (see §3.2). As can be seen from Figure 4 (upper left panels), typical relative 
errors of the broad components of H/3 and Ha are 6% and 3%, respectively. By relaxing 
the FWHM cut-off in H/? or Ha up to 2 200 km s~^, we could cover most of the bona fide 
NLSls with low spectral quality, at the expense of introducing a few spurious ones though. 
Imposing the stricter criterion of FWHM=2 000 km s'^ would reduce the sample size to 
1 885 objects. (2) The restriction of [OIII]/ H (3 < 3 was originally introduced to separate 
NLSls from Seyfert 2s where the resolution and/or S/N of the observed spectra were too low 
to do so. This is not the case for at least the majority of objects in our study, however, as 
discussed in §2.1. Consequently, this criterion becomes no long necessary. In fact, there is 
only one object (SDSS J104755. 93-1-073951. 2, see Figure 5 upper panels) that fails to meet 
this restriction. 

As a serendipitous finding, one object (SDSS J113323. 97-1-550415. 8) having the promi- 
nent, narrow "broad" Balmer emission fines [FWHM{H(3) ~ 2 100 km s"^ and FWHM{Ha) ~ 
1 900 km s~^\ characteristic of NLSl was actually a type 2 supernova. Its spectrum and 
SDSS image are shown in Figure 5 (lower panel). The object is located at the edge of the 
nearby blue compact galaxy UGCA 239 (SDSS J113323. 46+550420. 6) and has the same red- 
shift as the galaxy. It was noticed because its brightness decreased by ~ 3 magnitudes within 
less than one-year between the photometric and spectroscopic observations. Its disappear- 
ance from our later imaging observation of UGCA 239 using the 2.16m telescope of Beijing 
Observatory in January, 2005 confirmed its supernova nature. Accordingly we removed this 
object from our NLSl sample. 

3.2. Reliability of the Sample Parameters and Comparison with Previous 

Samples 

Our NLSls were drawn from a well-defined parent population of broad emission-line 
AGNs comprising candidates of both QSOs (Richards et al. 2002) and galaxies (Strauss 
et al. 2002 and Eisenstein et al. 2001). Considering the fact that the vast majority of the 
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sample were selected based on the optical selection criteria, we regard the sample simply as an 
optically selected sample. The emission-line parameters were measured with high accuracy 
by precise subtraction of the stellar and nuclear continuum and proper decomposition of the 
permitted emission-lines into a narrow and a broad component. Since the broad components 
of Ha and 11(3 were measured separately, we could check our emission-line fitting algorithm 
by comparing the two measurements. Figure 4 which plots FWHM(Ha) against FWHM(H/5) 
confirms the previously known correlation between the two. Our hnear fit, FWHM(Ha) = 
(0.842 ± 0.016) FWHM{HP), is consistent with the FWHM{Ha) = 0.873 FWHM{HP) 
found by Eracleous & Halpern (1994), and the 1 a deviation about 10% is only slightly larger 
than typical measuring errors (see the inset panel). We also compared the fluxes of the broad 
components of H« and H/9 and found F{Ha) = (3.028 ± 0.017) F{H 3) and a dispersion 
of ~ 0.36 in the Ha to H/? ratio (Figure 6). The good agreement between the independent 
measurements of the broad Ha and H/? components indicates that our emission-line fitting 
algorithm should be accurate and reliable. 

The method for identifying NLSls we used in this work is largely different from that 
used in WPM02; moreover, the WPM02 sample was drawn merely from the "QSO" database 
of the SDSS early data release (EDR). It would therefore be interesting to compare our 
sample with that of WPM02. We found that the 150 NLSl candidates in the WPM02 
sample are all included in our parent sample of broad emission-line AGNs. master sample. 
Our measured FWHM(H/3) are compared with those in WPM02 in the lower left panel of 
Figure 4. For those objects in which the narrow component of H/3 is neghgible, we found 
excellent agreement between the two measurements. However, for about 1/3 of the objects, 
our measurements of FWHM(H/3) are significantly larger than theirs. Going through these 
objects by visual inspection of their spectra revealed that the discrepancy is due to the non- 
negligible contribution of the H/? narrow component, which was not taken into account in 
WPM02. To demonstrate this effect we plot in Figure 4 (lower right panel) the FWHM(H/9) 
ratio of our measurements to theirs against the height ratio of the narrow to broad component 
of H/3, taken as an indicator of the relative strength of the narrow component. It can be 
seen that the measured FWHM(H/3) values in WPM02 start to deviate from ours when 
the narrow component becomes non-negligible, that the deviation gets worse as the narrow 
component increases, to reach a factor of 2 or more when the narrow component surpasses 
the broad component in height. 

As a result, 5 objects in the WPM02 sample were excluded from our NLSl sam- 
ple because their FWHMs of the H/3 broad component actually exceeded our criterion of 
2 200 km s~^. Their de-blended components of H/3 and Ha resulting from our analysis 
are shown in Figure 7. On the other hand, our sample includes 86 objects from the EDR 
database that were not included in the WPM02 sample. Of these, only 17 fall within the 
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range defined by the differing cutoffs, 2 000 km s''^ < FWHM(HP) < 2 200 km s-\ wfiile 
most of tliese because tliey were assigned by WPM02 to a separate category, Seyfert 1.5, by 
virtue of tlieir prominent narrow components of tlie Balmer emission-lines. Wc included these 
simply because their broad components had FWHM(Ha) or FWHM(H/3) that satisfied our 
cut-off criterion. We argue that the nature of these objects (see Figure 8 for some examples) 
should not be much different from classic NLSls, and there seem to be no particular reasons 
to exclude them based merely on the prominence of their narrow emission-lines. In addition, 
thanks to our starlight-nucleus decomposition and automated spectral fitting algorithm, our 
sample includes 13 NLSls drawn from the EDR Galaxy database, which was not considered 
in WPM02 at all; in fact, a total of 138 NLSls were culled from the Galaxy catalog from 
the SDSS DR3 database. The rest few objects not included in the WPM02 sample have 
either low spectral quality or large imccrtaintics in the measured spectral parameters. We 
accordingly believe that our NLSl sample has a greater degree of completeness than the 
WPM02 sample. 

The following caveats should be noted: Given the quality S/N 10 typical of the SDSS 
spectra for low redshift broad-line AGNs (BLAGNs), objects with very narrow "broad" 
permitted emission lines could not be effectively identified unless the Fell multiplets were 
sufficiently strong. Also, NLSls with large extinction may have been missed out. Observa- 
tions in multi- wavebands, such as the hard X-ray band, are needed to identify such objects. 

4. Sample Properties 

4.1. The Fraction of NLSls 

Our sample of NLSls is more than ten times the size of WPM02, while equally or even 
more complete and uniform. In this section we present some basic properties of the NLSl 
population, taking advantage of this unprecedented sample size. The NLSls as a whole make 
up 14% of all the broad emission-line AGNs. This fraction is consistent with the often quoted 
value of 15%, first suggested by Osterbrock (1988) and later confirmed by WPM02. Figure 9 
shows the distribution of NLSls and "normal" BLAGNs in redshift and nuclear luminosity 
(as represented by the g'-band absolute psf magnitude), as well as the NLSl fraction within 
overall BLAGNs. The NLSl fraction appears to increase slightly up to a redshift of 0.45, 
after which it decreases dramatically toward high redshift. This is actually not a reflection of 
the cosmic evolution of the NLSl fraction but is rather one of its dependence on the nuclear 
luminosity (as seen in Figure 9), which is strongly correlated with redshift for a flux- limited 
sample. 
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The fast decrease of the NLSl fraction with increasing nuclear luminosity in the high 
luminosity range results actually from the emission line width-luminosity relation, which is 
discussed below in § 4.2. The strong positive correlation of the NLSl fraction with luminosity 
in the low luminosity range Mg < —22™ is unexpected. The fraction peaks at Mg ~ —22™, 
where NLSls make up about 20% of all SDSS BLAGNs; this peak apparently results from 
a dip in the luminosity distribution of BLAGNs here (see Figure 9), without there being a 
corresponding one in that of NLSls. It is hkely that this dip in the BLAGNs distribution 
is due to selection effects, in the sense that most of low luminosity AGNs were targeted 
as "galaxies" while high luminosity ones as "QSOs" in the SDSS. However, we argue that 
the peak in the NLSl fraction should be real because objects of both the NLSl and BLSl 
types were selected according to the same criteria in the SDSS. It may be a result of the 
combination of the local black hole mass function and the distribution of the mass accretion 
rate. Detailed discussion on this issue is beyond the scope of this paper and will be presented 
elsewhere. 

It was once thought that NLSls arc radio quiet; however, the past few years saw the 
discovery of about a dozen "radio-loud" NLSls (Siebert et al. 1999; Grupe et al. 2000; Zhou 
& Wang 2002; Zhou et al. 2003). Two objects among these, both identified from the SDSS, 
are very-radio-loud (Zhou et al. 2003; Zhou et al. 2005) according to the criterion given 
in Sulentic et al. (2003). The incidence of radio- loud NLSls, or the dependence of the 
fraction of NLSls on the radio properties is of particular interest. This is because radio- 
loud NLSls show some extreme characteristics opposite to typical radio-loud AGNs, e.g. 
the extremely steep soft X-ray spectrum in contrast to the usual flat spectrum (Sulentic et 
al. 2000). This issue can be examined in more detail with our large NLSl sample. Out 
of the 2 Oil NLSls, 142 were detected in the FIRST (the Faint Images of the Radio Sky 
at 20 cm survey, Becker et al. 1995). The radio detection rate is 7.1 ± 0.2%, which is 
significantly less than what we found for BLSls at 10.2 ± 0.1%. It is remarkable that all the 
radio sources of the NLSls detected in the FIRST were unresolved in their radio images. 
The distribution of radio-loudness R, (defined as the logarithm of the flux ratio between 1.4 
GHz and the optical g-band) is shown in Figure 12 (left panels) respectively for the NLSls 
and BLSls detected in the FIRST. Compared to the BLSls, the radio-loudness distribution 
of the NLSls drops much faster beyond i? ~ 1, indicating a decline of the NLSl fraction 
toward high R values. The NLSl fraction is ~ 15.5% in radio-quiet and radio-intermediate 
AGNs (i? ^ 1), similar to that from other optically selected samples; it falls down to ~ 10% 
in moderately strong radio AGNs (1 ^ -R ^ 2) and further down to ~ 5% in powerful 
radio AGNs (i? > 2). In particular, we found only 5 very radio-loud NLSls in the whole 
sample, including the two we reported previously, SDSS J084957.98+510829.0 and SDSS 
J094857. 32-1-002225. 5. These two very radio-loud NLSls show flat spectra and dramatic 
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variability in the radio band, indicating the presence of relativistic jets pointing toward the 
observer. In fact, SDSS J084957.98+510829.0 had been observed before the SDSS and was 
classified as a BL Lac object with a wrong redshift (Zhou et al. 2005 and references therein). 
For the other three sources no significant radio variabihty was found in time spans of several 
years between the NVSS (the NRAO VIA Sky Survey, Condon et al. 1998) and the FIRST. 
With a spectral index of — 0.4, the radio spectrum of the most powerful NLSl, SDSS 
J104732. 68+472532.1, is somewhat steeper for a blazar-hke source. The dependence of the 
NLSl fraction on radio power (Figure 12, right panels) is similar to that on the radio-loudness. 
A detailed account on the radio properties of NLSls will be presented in a forthcoming paper 
of the series. 

Historically, X-ray surveys have been an important tool for the finding and study of 
NLSls. For instance, Puchnarewicz et al. (1992) found that about 50% of their soft X-ray 
selected AGNs were NLSls; a similar result was also obtained by Grupe (2004). However, 
we found that the overall NLSl fraction is only ~ 19.4 ± 0.8% for our optically selected 
sample. Plotted in Figure 13 (left panel) is the dependence of the NLSl fraction on the 
soft X-ray fiux density for the BLAGNs with X-ray counterparts in the RASS; it reveals a 
strong dependence — the NLSl fraction increases dramatically with the X-ray fiux. At the 
faint end, the NLSl fraction is ~ 15%, similar to that found in optically-selected BLAGN 
samples. At the bright end, it increases up to ~ 40%, close to that of the soft X-ray selected 
AGN samples reported previously; this is where the NLSl fractions for optically selected and 
soft X-ray selected AGN are reconciled. This is not surprising, as above this high flux level 
optically selected and X-ray selected BLAGN samples are largely overlapping. Interestingly, 
we also found a trend of decreasing fraction of NLSls toward high X-ray luminosity (Figure 
13, right panel). 



4.2. The Luminosity— Emission Line Width Relation 

It has been claimed by several authors (e.g.. Miller et al. 1992; WGOl) that the nuclear 
luminosity is correlated with the width of the broad B aimer emission-lines for BLAGNs, 
including NLSls. We plot in Figure 10 line luminosity versus line width for Ha and H/9, 
respectively, for our sample. A weak yet significant correlation is present between the two 
parameters, despite the small range of the line width (the Spearman coefficient r = 0.29 for 
FWHM(H/?) and L(H/?), and r = 0.34 for FWHM(Ha) and L(Ha), corresponding to chance 
probabilities much less than 0.1% in both cases). The correlations appear to result from the 
presence of an upper boundary, as indicated by the dashed lines in Figure 10. As will be 
seen below, the limit is likely caused by the existence of upper limits in the accretion rate 
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in units of the Eddington rate. 

The central black hole mass can be inferred from the size and the velocity dispersion of 
the broad emission-hne region (BLR), 

Mbh oc Rblr V\ (2) 

provided the BLR is viriahzed. We can use the width (FWHM) of Balmer emission-lines to 
estimate the velocity dispersion of the BLR. Supposing all AGNs have on average the same 
ionization parameters, the same ionizing spectral energy distribution (SED), the same BLR 
densities and column densities, then we would expect a simple scaling relation between the 
BLR size and the monochromatic luminosity at 5100 A (e.g., Peterson 1993), 

Rblr « [ALa(5100)]°-^ (3) 

Hence we can express the virial mass of the central black hole mass as 

Mbh oc [ALa(5100)]°-^ [FWHM{Hp or Ha)^. (4) 

For our NLSls, we find that AL;^(5100) is tightly correlated with the luminosity of H/3 broad 
component, L(HP) (Figure 11, left panel). A linear fit yields, 

%{ALa(5100)} = (8.57 ± 0.28) + (0.840 ± 0.007) log{L{H(3)}. (5) 

Also, assuming that all objects have the same SED, we get 

Lboi oc ALa(5100) oc [L{Hp or Ha)]^-^^^. (6) 

Therefore, for objects emitting at the Eddington luminosity, 

Lboi = LEdd oc Mbh, (7) 

a given value of FWHM (Hp or Ha) corresponds to a value of L{HP or Ha), 

L{H(3 or Ha) = C [FWHM{H(3 or Ha)]^-'^^^. (8) 

The exact value of the constant C depends on the shape of the ionization continuum, the 
geometry and physical condition of the BLR, and the cosmology adopted. This correspon- 
dence is visuahzed in Figure 10 by the dashed hne with the constant C so chosen that L{,o; 
is a few units of LEdd- It can be seen that most of our NLSls populate below the line. This 
result implies that the radiation of the majority of AGNs including NLSls cannot be much 
higher than a few units of the Eddington luminosity. The correlation between the broad 
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line width and the nuclear luminosity naturally explains the decrease of the NLSl fraction 
toward high luminosities reported in §4.1. 

A strong correlation between the luminosity of continuum and [OIII] emission line was 
also found for our NLSls (see right panel of Figure 11), though not as tight as the L{HP) — 
ALa(5100) correlation. The best linear fit to the data points yields, 

/o^{ALa(5100)} = (7.11 ± 0.56) + (0.885 ± 0.0133) logmoiii]}- (9) 

This relation is useful for estimating the intrinsic luminosity of narrow emission-line AGNs 
based on [OIII] luminosity (e.g., Zakamska et al. 2003). If the central black hole mass of 
these type 2 AGNs can be reliably estimated, one can evaluate their mass accretion rate and 
search for type 2 counterparts of NLSls in those objects with high Eddington ratios L/ Lem- 

4.3. Broad Emission Lines 

Though not a defining property, NLSls usually show strong optical Fell emission fines. 
Out of the 2 Oil NLSls, the Fell multiplets were detected in 1 787 objects at the > 3 cr 
level. For the rest the SDSS spectra were too noisy to yield reliable measurements and the 
3 a upper limits were calculated (see Table 1). The relative strength of the Fell multiplets 
is conventionally expressed as the flux ratio of Fell to H/?: /?,i57o = p<^^^^^^-^^-^*->^'^ ^ where 
Fe//AA4434 — 4684 denotes the flux of the Fell multiplets integrated over the wavelength 
range of 4434 - 4684 1, and H/3 the flux of the broad component of H/3 (e.g., VVGOl). The 
distribution of -R4570 for the NLSl sample is shown in Figure 14. The average is -R4570 = 0.82 
and the dispersion at 1 cr is 0.37 when only the 1 787 objects with reliable Fell detection were 
considered. This is significantly larger than the typical value of i?457o ~ 0.4 found in normal 
AGNs (Bergeron & Kunth 1984). The inclusion of the non-detections would not change the 
results significantly, since they only account for ~ 10% of the whole sample. About a quarter 
of all the NLSls are moderately strong Fell emitters with -R4570 ^ 1, while this fraction is 
only ~ 5% in normal AGNs (Lawrence et al. 1988). We identified 17 new extreme Fell 
emitters (i?457o > 2) from the present NLSl sample, which doubled the number of known 
objects of this kind (VVGOl; Zhou et al. 2002; Zheng et al. 2002). 

The merits of our sample are its large size and uniformity, thanks to which some weak 
correlations are expected to be identifiable. Hence, the sample is suited for an examination of 
whether the previously known correlations for broad fine AGNs can be extended to smaller 
line widths. The El space involves primarily the correlations between three parameters: 
FWHM(H/3), -R4570, and soft X-ray spectral index, F. In Figure 15 -R4570 is plotted against the 
FWHM of the broad components of H/3 and Ha, respectively. We found that the known anti- 
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correlation between -R4570 and the FWHM of the broad Balmer hne, though weak, continues 
to exist in NLSls. The Spearman correlation coefficient is r=-0.23 for -R4570 vs FWHM(H/3) 
and r=-0.31 for -R4570 vs FWHM(Hq;), corresponding to chance probabihties of <^ 0.1% 
in both cases. The correlation is slightly stronger between -R4570 and FWHM(Ha) than 
between -R4570 and FWHM(H/3); this is because, in general, the measuring uncertainties are 
smaller in FWHM(Ha) than in H/5. GaskeU (1985) suggested that the large -R4570 in NLSls 
was due to their Balmer emission lines being weak rather than their Fell emission being 
strong; this was confirmed by VVGOl by showing a moderately strong correlation between 
H/3 luminosity and -R4570 found in some 200 broad line AGNs including 64 NLSls. We did 
not find such a correlation for NLSls alone, however. This result indicates that NLSls form 
a distinct subclass, at least as regarding those showing rather broad low-ionization lines, 
namely "population B" (Sulentic et al. 2000). 

Interestingly, strong correlations were found between the luminosity of the Balmer emis- 
sion lines and the equivalent widths of both the Balmer and Fell emission lines (the "inverse" 
Baldwin effect) for NLSls. We do not know exactly the origin of this strong correlation, which 
is rather weak — if exists at all — in BLSls. Contamination from stellar light can be ruled out 
as we believe that starlight has been properly subtracted out. Either a luminosity-dependent 
SED, for which the optical to UV spectrum flattens as luminosity increases, or an increase 
in the covering factor of BLR with increasing luminosity may explain such a correlation. 

A strong correlation between the ROSAT photon index F and FWHM(H/3) has been well 
established for broad emission-line AGNs (Laor et al. 1997; Wang, Brinkmann, & Bergeron 
1996; Boiler et al. 1996). While BLSls show a rather small scatter in their photon indices, 
NLSls show a larger dispersion of F ~ 1 — 5 (Zhou & Wang 2002). It is not yet clear 
whether the r-FWHM(H/3) correlation could extend to very small Balmer line widths. A 
total of 635 NLSls have X-ray counterparts in the RASS catalog; however, only a fraction 
have enough photon counts for an estimation of the photon index. We searched for more 
data from the ROSAT Source Catalog of pointed observations (RSC, Voges et al. 1999) in an 
attempt to improve the X-ray data quality. X-ray sources that lie within 30" of the optical 
positions were regarded as real matches. In all 136 RSC sources were found (all detected in 
the RASS) and their RSC data were used in preference to their RASS data, because of their 
longer exposures. Following Schartel et al. (1996) and Yuan (1998), we calculated the X-ray 
photon indices and X-ray fluxes using the two hardness ratios and the count-rates given in 
the ROSAT catalogs for an assumed Galactic absorption column density. 

Figure 17 plots the ROSAT photon indices F versus the FWHM of H/3 and Ha for the 310 
NLSls, for which F could be estimated with reasonable certainty. Of particular interest, there 
appears to exist a turn-over in the correlation trend at FWHM around 1 000 km for both 
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H/3 and Ha. For FWHMs above the turn-over (FWHM > 1 000 km s"'^) T is significantly 
anti-correlated with both FWHM(Ha) and FWHM(H/?); applying the Spearman test gave 
a correlation coefficient r = —0.24 and r = — 0.2f , and a chance probability of p = 0.1% 
and p = 0.04%, for r-FWHM(Ha) and r-FWHM(H/3), respectively In the FWHM range 
below 1 000 km on the contrary the trend is reversed, with r = 0.366, p = 1.7% 
and r = 0.316, p = 1.2% for the r-FWHM(Ha) and r-FWHM(H/3) relations, respectively 
Objects with the steepest spectral slopes (the peak in the F distribution) were found to 
have FWHM values around ~ 1 000 km s'^. For the 61 NLSls in the FWHM range of 
800 - 1 200 km s'^, the mean photon index is < F >= 2.9 ± 0.4, with 43 objects (~ 70%) 
have F > 2.75 and only 2 objects have F < 2.0. 



4.4. Nctrrow Emission Lines 

For 751 NLSls in our sample the fluxes of the H/3 narrow component could be rehably 
measured at the > 3 cr level; for most of the rest meaningful upper limits could be obtained. 
We plot in Figure 18 (upper panel) the distribution of the flux ratio of [OIII]A5007 to 
the narrow component of H/3, [OI I I]X5007 / H P'^ . For those objects with only upper limit 
on the H/9 component, a lower limit of [0 1 1 1]X5007 / H f3^ is set at 10 . Objects having 
[0 1 1 1]X5007 / H (3^ < 5 were found in only ~ 15% of the NLSls of our sample; and most of 
those NLSls with low [OIII]\5007/Hf3'^ ratios are actually "composite" objects in which a 
large fraction of the emission lines come from HII regions in the host galaxies. This result 
is consistent with that of VVGOl but is contrary to that of Rodriguez- Ardila et al. (2000). 
We plot in Figure 18 (lower panel) the diagnostic diagram for the NLR of 168 NLSls, for 
which the two sets of narrow emission hues, [OIII]/ H 13'^'^ and [Nil]/ Ha^"^ could be rehably 
measured. About 10% of these NLSls are scattered into the region occupied by star-forming 
galaxies; this fraction is also consistent with that of VVGOl (6/64 in their sample). One 
important implication is that the majority of the type 2 counterparts of NLSls are actually 
Seyfert 2 galaxies. 

Using a sample of ~ 400 broad emission-line AGNs selected from the first data release 
of the SDSS, Boroson (2005) found that the low-ionization forbidden lines ([Oil], [Nil] and 
[SII] doublets) had mutually consistent redshifts, while in about half of the AGNs the [OIII] 
doublet was blueshifted with respect to the redshift defined by the low-ionization forbidden 
lines. The average blueshift of [OIII] was 40 km and the largest shift could reach 
> 400 km s~^. This phenomenon also appears in our NLSl sample (Figure 19, panels a 
and b). The peaks of the [OIII] hue in ~ 60 NLSls are shifted by > 400 km s'^. The 
spectrum of one such object is shown in Figure 20 (upper panel). These so called "blue 
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outliers" (Zamanov et al. 2002) can provide a valuable tool in the study of the dynamics of 
NLR and its possible connection with other AGN components. 

We also found a few NLSls showing double-peaked profiles of narrow emission-hnes 
(Figure 20, lower panel). Possible explanations of the origin of such profiles invoke bipolar 
outflows, a disk-hke emission hne region, and binary black holes with separated NLRs (Zhou 
et al. 2004). 

4.5. Notes on Other Rare Type NLSls 

One remarkable property of this large sample is its diversity. Rarities could be readily 
identifled, such as very radio-loud NLSls, "blue outliers", and objects with double-peaked 
narrow emission lines as mentioned above. Below we present several other rare types of 
objects, namely "UV-deficient" NLSls, "Fell-deflcient" NLSls, and NLSls with broad ab- 
sorption line troughs (BALs). As illustrative examples, some of their spectra are shown in 
Figure 21 and 22 and 23. Rare type objects are very useful for examining AGN models. 

It is generally believed that each NLSl harbors a relatively small black hole accreting 
at a very high accretion rate. According to this scenario, NLSls would have a flat spectral 
shape in the optical-to-UV band. While this was found to be indeed the case for the majority 
of the sample, some NLSls show deficient radiation in the blue part of their SDSS spectra. 
Most of these "UV-deficient" spectra can be explained as due to heavy extinction caused by 
a presumed dusty torus and/or dust lane in the host galaxies. In a few objects, however, 
this unusual spectrum cannot be explained by extinction. For example, the UV dip in the 
spectrum of SDSS J142033. 7-1-573900. 9 (Figure 21, upper panel) cannot be entirely accounted 
for by extinction because the flux ratio of H^/HP — 0.40 ± 0.02 is close to the theoretical 
value expected for a spectrum free from extinction. 

Another rare class of NLSls is concerned with the strength of the Fell multiplets. 
While the Fell emission lines are generally found to be strong in most NLSls, they are very 
weak, or even undetectable, in some objects with very narrow "broad" Balmer emission lines 
(Figure 22). The closest type to these objects is dwarf Seyfert 1 galaxies, such as NGC 4395 
(e.g., Filippenko & Sargent 1989). The difference between the two classes lies in that the 
Fell deficient NLSls found here have a relatively high luminosity with respect to the black 
hole mass, and hence a very high accretion rate. 

NLSls and broad absorption line (BAL) QSOs share many common observational prop- 
erties (see Brandt & Gallagher 2000 for a review). This similarity motivated several authors 
to investigate possible connections between the two. For instance, there is the suggestion 
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that BAL QSOs are misaligned NLSls (Sulentic et al. 2000). We have identified 6 NLSls 
with possible BAL (3 of the candidates were serendipitously discovered in the SDSS DR4), 
whose SDSS spectra are shown in Figure 23. The facts are perplexing. On one hand, BAL 
QSOs are thought to be otherwise normal QSOs seen nearly edge-on according to the pop- 
ular unification models. On the other hand, there is convincing observational evidence for 
a small inclination angle for at least some NLSls (e.g., Zhou et al. 2003). Despite many 
observed properties of BAL QSOs can be explained by orientation-based unification models, 
evidence begins to mount that polar outflows do exist, at least in some BAL QSOs (Zhou et 
al. 2006). It appears that the real situation is much complex than is assumed by the simple 
models. 



5. Black Hole-Bulge Relation in NLSls 

5.1. Mbh — cr* Relation 

A growing body of evidence accumulated in the past decade suggests that the evolution 
of supermassive black holes (SMBHs) and that of their host galaxies are tightly linked. Ko- 
rmendy & Richstone (1995) found that SMBH mass, Mbh, is correlated with the mass of 
the spheroidal component, Mhuige, in seven nearby galaxies. This relationship was later con- 
firmed and quantified by Magorrian et al. (1998) based on a sample of 36 nearby galaxies with 
Hubble Space Telescope (HST) photometry and ground-based kinematics data. Ferrarese & 
Merritt (2000) and Gebhardt et al. (2000) demonstrated that Mbh is tightly correlated with 
the stellar velocity dispersion, cr*. Nelson et al. (2004) measured the bulge stellar velocity 
dispersion in 14 Seyfert 1 galaxies whose SMBH masses were determined using the reverber- 
ation mapping technique, and showed that the Seyfert galaxies followed the same Mbh — cr* 
relation as non-active galaxies. The tight correlations for both normal and active galaxies 
indicate that the most of the SMBH mass was built up in the past by accretion of the same 
gas that also formed the stars in the host galaxy spheroid. It is of particular interest as to 
when and how SMBH and the surrounding spheroid of stars settled down at the right masses. 
A straightforward approach to investigate the physical link between the SMBH growth and 
the bulge formation is to measure the SMBH masses and bulge properties in high redshift 
AGNs, and compare them with those in low redshift AGNs and quiescent galaxies. This is a 
very hard undertaking given the observational difficulties. One alternative is to study young 
AGNs in the local universe. It has been suggested that NLSls are normal Seyfert galaxies 
in their early stage of evolution, when the SMBH is growing fast via accreting mass at a 
high rate (e.g., Wang, Brinkmann, & Bergeron 1996; Boiler, Brandt & Fink 1996; Laor et al. 
1997, Mathur 2000). In this sense NLSls may serve as a surrogate for high redshift AGNs 
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(Wang & Lu 2001). 

A direct probe of the Mbh — o'* relation is practically useful since, in principle, the stellar 
velocity dispersion of a bulge in nearby active galaxies can be determined more accurately 
than the bulge luminosity or compactness (Haring & Rix 2004). In practice, however, width 
of the [OIII] emission line is often used in stead of stellar velocity dispersion due to observa- 
tional difficulty of the latter. Using such an approach several groups of authors explored the 
Mbh — <7* relation for NLSls recently, with conflicting results: Mathur et al. (2001), Grupe 
& Mathur (2004), Bian & Zhao (2004), and Botte et al. (2004) found that NLSls showed 
a systematically lower Mbh /Mbuige ratio than normal broad line AGNs, while Wang & Lu 
(2001), Wandel (2002; 2004) did not find a clear difference between the two. The difference 
may be attributed, at least partly, to the different line widths used by the different authors. 
Wang & Lu (2001) argued that only symmetric line component can trace the gravitational 
potential of the bulge and so they used the width of the blue- wing subtracted [OIII] line from 
VVGOl; in contrast, most others adopted the total [OIII] line profile. It has been shown 
that the [QUI] line profiles in most Seyfert galaxies are evidently asymmetric, so signifying 
complex kinetics of the emission line regions (VVGOl and references therein). Moreover, 
contamination from the Fell emission aggravates the problem, which is especially true for 
NLSls. To remove all these uncertainties from the study of the SMBH-bulge relation, it is 
essential to have direct measurements of cr* in stead of the [OIII] line width. The underly- 
ing physical justification is clear, since the stellar kinematics is more indicative of the bulge 
gravitational potential than is the gas in the NLR of the AGN. As an attempt along this line, 
Botte et al. (2005) obtained the near-infrared spectra of 8 NLSls with a moderately high 
resolution of ~ 50 km per pixel, and measured directly the stellar velocity dispersions 
using the CallAA 8498,8542,8662 absorption triplet. They found that (J[oiii] ^ ioi 8 of 
the 10 NLSls with cr* available so far, and demonstrated that the NLSls follow the same 
Mbh ^ o"* relation as normal broad line AGNs and non-active galaxies. However, the sample 
of Botte et al. (2005) is too small to give conclusive results, as was warned by the authors 
themselves. A large NLSl sample with properly measured u* is utterly needed. 

As discussed above in §2.2, in the process of decomposing the stellar and nuclear com- 
ponents in NLSls, cr* can be obtained automatically^. Though the resolution and S/N ratio 
of the SDSS spectra are not perfect for exploring the Mbh — cr* relation in NLSls, the 
present sample is so large as to render it sensitive for identifying any weak trend. Apart 
from the spectral resolution and signal-to- noise ratio, the measurement uncertainty of cr* 
is also strongly dependent on the relative contribution of the stellar light. We found that 



^Because the spectral resolution of the galaxy templates we used is similar to that of the SDSS, decon- 
volution of the SDSS spectra was not necessary to obtain the stellar velocity dispersions. 
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(7* could be measured with acceptable errors for spectra with a median signal-to-noise ra- 
tio S/N > 10 and a relative contribution of starlight more than 40% of the total flux at 
5100A. A small number of objects have measured a^, values as small as < 50 km 
which is actually beyond the resolution of the SDSS; they were discarded from the analysis. 
Adopting these criteria yielded 308 objects with, we believe, reasonably well determined a^,; 
they formed the subsample which we used to test the black hole-bulge relation for NLSls. 
Figure 24 plots, for this subsample, the estimated Mbh versus cr*, with the Mbh determined 
using the hne width-luminosity mass scaling relation (e.g., Kaspi et al. 2000; McLure & 
Jarvis 2002; Dietrich & Hamann 2004), 

Mq ^ 10^4 erg s'^^ ^ 1 000 km s'^ ^ ^ ^ 

The well known Mbh^(^* relation for normal galaxies, as parameterized by Tremaine et al. 
(2002), 

/OS(^) = (8.13± 0.06) + (4.02 ± 0.32) (H) 

is marked in the plot by the dashed line. A correlation is clearly present, which follows 
the trend of the Mbh^<^* relation for galaxies (the Spearman correlation coefficient r=0.27 
giving a chance probability of p = 2 10~^). Meanwhile, the data points of NLSls show 
a systematic downward shift relative to the relation for galaxies as defined in Equation 5, 
with the bulk of the subsample (~ 90%) located below the dashed line. The results indicate 
that, firstly, the correlation between Mbh and a* remains valid for NLSls; and secondly, the 
SMBHs in NLSls are underage and their growth lags behind the formation of the bulges, as 
suggested by e.g., Mathur et al. (2001). We stress that our results are actually consistent 
with Botte et al. (2005); their much smaller sample prevented the authors from drawing 
the same conclusion. In fact, only 2 of the 10 NLSls in their sample are located above the 
Mbh-c* line for normal galaxies and normal AGNs. 

Caution should be exercised when accepting the above results, considering that the 
derived cr* may be affected by contamination by the rotation of the galactic disk. This 
issue presents a matter of concern here since, for NLSls at typical redshifts in our sample, 
the fiber aperture would cover the whole galactic disk. For a further examination of the 
reality of the Mbh — o"* relation obtained above, we tried to minimize the contamination 
effect of the galactic disk rotation. We visually inspected the SDSS images of NLSls at 
redshifts < 0.1 with directly measured cr*, and picked out objects for which either the host 
galaxies appear to be face-on or the 3" fiber aperture is dominated by the galactic bulge 
contribution. This selection yielded 33 objects; their SDSS images are shown in Figure 25. 
For these objects, we believe that the contamination of the disk, even if significant, should 
not lead to overestimation of cr*. Their Mbh — cr* relation is plotted in Figure 26. It can be 
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seen that the same trend of correlation persists, and all the data points are distributed well 
below the expected Mbh — relation for galaxies of Tremaine et al. (2002). 

5.2. cr*— Narrow Line Width Relation 

Since we have relatively good measurements of cr* for the subsample, we could use 
them as a standard calibrator to test potential correlations with the width of any narrow 
emission-lines. This has practical potential in the sense that in many NLSls, especially 
luminous ones, the measurement of a* is not available, that maybe the line width could 
serve as a good surrogate. Below we test the previously reported relations between cr* and 
the FWHM of the [Nil] and [OIII] lines. We plot in Figure 27 the FWHM[NII] values 
(deconvolved to remove the effect of the SDSS spectral resolution) versus cr* for 240 out of 
the 308 NLSls, for which the [Nil] doublet is detected at the > 10 cr level. A significant 
correlation was found at a probability level p = 1.8 10~^^ (the Pearson correlation test). The 
objects were found to be distributed almost symmetrically with respect to the often quoted 
relation FWHM[NII] ~ 2.35 (7* with a la deviation of ~ 30% (solid hue in Figure 27). The 
result indicates that the [Nil] emission hue can be used as a surrogate for stellar velocity 
dispersion in NLSls. We extended the above test to a larger sample of Seyfert 2 galaxies 
compiled from the SDSS, comprising some 3 000 objects with high spectral S/N ratios. A 
correlation is found very similar to that for the NLSls (Figure 28). A linear fit yielded 
FWHM[NII] ~ 2.62 a,, as compared to the often quoted FWHM[NII] ~ 2.35 a,. We 
conclude that, for Seyfert 2s, FWHM[NII] can also be used as a surrogate of cr*. We find 
that, after taking all the scatter into account, FWHM[NII] can predict black hole masses 
to within a factor of 3. 

In comparison, the [OIII]-cr* correlation is only marginal. Our result suggests that the 
[OIII] line width, which has been widely used in the literature, is actually not a good indicator 
of the (7* and, hence, of the black hole mass. In fact, this is not unexpected considering that, 
in more than half of the NLSls in our sample, the [OIII] lines have complex profiles. In a 
systematic investigation of the stellar and gaseous kinematics in narrow emission-line AGNs, 
Greene & Ho (2005) also noted the different behavior of the low-ionization lines and [OIII], 
and pointed out that cr^oiii] cannot be directly used to replace o"*. Using a two-component 
Gaussian fit to the [OIII] line profiles, the authors found that the core of the [OIII] line could 
trace a^, statistically but with a considerable scatter. This approach, however, cannot yield 
meaningful results for most of the objects in our sample because of the low S/N ratio and/or 
the effect of the residuals of the Fell emission in their SDSS spectra. 

Having established the good FWHM[NII]-a^. relation, we can extend the above ex- 
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amination to the black hole-bulge relation by making use of FWHM[NII]. The advantage of 
so doing is that, compared to o"^,, FWHM[NII] is much easier to measure and is available for 
a larger number of NLSls. In Figure 29, we plot Mbh versus FWHM[NII] for 613 NLSls 
with reliably measured FWHM[N11]. A correlation is present, but there is a systematic de- 
viation from the equivalent relation expected for galaxies given by Tremaine et al. (2002) 
transformed from FWHM[NII] = 2.35 cr* (dashed line). The results confirm those obtained 
above for a smaller sample by making direct use of (Figure 24). 



6. Conclusions and Prospects 

6.1. Summary of Conclusions 

We have systematically analyzed the spectroscopic data from the QSO and galaxy 
database in the SDSS DR3. After proper subtraction of the stellar and nuclear continua 
as well as the Fell multiplets, prominent emission lines were carefully modeled and mea- 
sured with typical relative errors less than 10%. A large sample of ~ 2 000 NLSls has been 
selected based on a simple well-defined criterion: the presence of a "broad" component of 
the H/? (or Ha) line with FWHM < 2 000 km s~^. The present sample outnumbers pre- 
vious NLSl samples by a factor of ~ 10. Taking advantage of such an unprecedented large 
and uniform sample with accurately measured spectral parameters, we carried out various 
statistical analyses, some of which were possible only for the first time. 

Overall, the often quoted frequency ~ 15% is confirmed for finding NLSls in optically 
selected broad emission-line AGNs; meanwhile, the frequency is also found to be strongly 
dependent on the optical luminosity. The fraction of NLSls peaks around Mg ~ —22"^ 
(the SDSS g'-band absolute magnitude) and drops quickly toward both the high- and low- 
luminosity ends. We interpret the lower chance of finding NLSls in high luminosity AGNs 
as a result of an imposed upper limit on the Eddington ratio (L/LEdd), of a few units. 
The smaller fraction in the low luminosity range may arise from the combination of two 
factors: the shape of the SMBH mass function and the distribution of the accretion rate. 
Low luminosity NLSls contain a very small SMBH mass accreting at a very high rate, which 
are rare in the local universe. We do not confirm the previous results that ~ 50% of soft 
X-ray (e.g. ROSAT) selected AGNs are NLSls. Such a fraction is found only for X-ray 
bright AGNs, but drops quickly with decreasing soft X-ray flux. The frequency of flnding 
NLSls in faint ROSAT AGNs (< IQ-^^.s ^-i ^^-2 ^j^g 0.1-2.4 keV band) is actually 
the same as in optically selected AGNs. The dependence of the NLSl fraction on the soft 
X-ray luminosity shows the same trend as on the optical luminosity, and the fraction peaks 
around ~ 10^^-^ erg s~^. Deeper X-ray surveys than the RASS are needed to pin down the 
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exact peak luminosity. Moreover, we find that the NLSl fraction is also strongly dependent 
on the radio- loudness and the radio luminosity. The frequency of finding NLSls in radio- 
quite AGNs (with logarithmic radio-to-optical flux ratio i? < 1) is ~ 15%, similar to that 
for optically selected AGNs. This value drops to ~ 10% in moderately strong radio AGNs 
(1 < i? < 2), and further down to ~ 5% in powerful radio AGNs (i? > 2). 

Wc find that on average the relative strength of the Fell emission (-R4570) in NLSls is 
about twice of that in normal AGNs; -R4570 is anti-correlated with the width of the broad 
component of the Balmer emission lines. The well known correlation between the width 
of the broad low-ionization line and the soft X-ray spectral slope is found to hold true for 
NLSls in general; in addition, it shows some quite peculiar behavior: F is anti-correlated 
with FWHM(H/5) in the range FWHM{H/3) > 1 000 km s'^, then the trend is reversed in 
the range FWHM{H(5) < 1 000 km s'^ . We find that the equivalent width of the H/3 and 
Fell emission lines are strongly correlated with the H/3 and the continuum luminosity. We 
do not find any difference in the NLR properties between NLSls and normal AGNs. About 
10-20% of NLSls show abnormally low line ratios of [OIII] to the narrow component of H/3, 
which is a result of the contribution of HII regions in the host galaxies. Blueshifts of the 
[OIII] emission line were found in about half of the NLSls, similar to that in normal AGNs. 

Stellar velocity dispersions o"* could be measured within statistically meaningful uncer- 
tainties for 308 NLSls. For this subsample we examined the Mbh^<^* relation for NLSls. 
We found that the black hole mass is also significantly correlated with the stellar velocity 
dispersion; however, the bulk of the NLSls have black hole masses statistically smaller than 
expected from the Mbh — o'* relation for normal galaxies and AGNs. This results support the 
hypothesis that NLSls are underage AGNs, and the growth of their SMBH lags behind the 
formation of the galactic bulge. With a typical relative error ~ 30%, the line width of [Nil] 
is found to be a better indicator of cr* than that of [OIII] . By using a larger subsample with 
reliably measured FWHM[NII], the Mgij-bulge relation for NLSls was further examined 
and similar results were obtained. 



6.2. Remarks on Future Work 

Limited by the scope of the this paper, most of the analyses presented above are only 
phenomenological. Systematic and detailed investigations based on this large data set of 
some of the interesting aspects as outlined in the Introduction (§ 1) are just underway. For 
instance, we have initiated a program to study the X-ray spectral and temporal properties 
of the NLSls in our sample by making use of archival serendipitous data from ROSAT, 
XMM-Newton and Chandra. The multiwavelength SED can be constructed and studied 
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by compiling data from the existing archives in various wavebands. The present results 
on the Mbh — c* relation for NLSls are only statistically meaningful; observations with 
higher spectral resolution and S/N ratios are needed for more precise results. Moreover, 
mounting evidence accumulated in the past decade indicates that the AGN phenomena and 
the star formation history are closely linked; the preliminary result on the Mbh — o"* relation 
for NLSls may represent an important link in the sequence of galaxy evolution. In this 
direction, we are going to study the stellar content of the host galaxies in the hope of finding 
any clues to the starburst-AGN connection. 

Some peculiar objects of various rare types have been identified in our NLSl sample, 
namely, very-radio-loud NLSls, [OIII] "blue outliers", double-peaked narrow line emitters, 
UV-deficient NLSls, and BAL NLSls. Detailed studies of their properties, which are ex- 
pected to be useful for constraining AGN models, will be addressed in forthcoming papers. 

In the above and other related studies of NLSl, a major source of uncertainty comes from 
the determination of the black hole masses. Mass estimation using a one-epoch spectrum is 
subject to large uncertainties (cf., Vestergaard 2004); reliable measurements of the SMBH 
mass are limited to only a few NLSls so far, however. With interest we noted that a few 
percent objects in our NLSl sample have duplicated observations in the SDSS, some of which 
show large amplitude variability (see Figure 30). This property presents us a great potential 
of measuring the black hole masses of NLSls by making use of the much more accurate 
reverberation mapping technique. We plan to monitor a NLSl sample with a moderate size 
to select candidates and then to carry out the reverberation mapping. The work will also 
enable us to explore the dynamics of the BLR for NLSl. 
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Table 1. Emission-Line Properties of SDSS NLSls (a portion for guidance of its form and content). Col. 1: object name in J2000, 
col. 2: redshift given by the SDSS spectroscopic pipeline, col. 3: monochromatic luminosity at 5100 A, col. 4: featureless nuclear 
continuum fraction at 5100 A, col. 5: H/3 narrow component flux, col. 6: H/? broad component flux, col. 7: HI3 broad component 
FWHM, col. 8: [OIII]A5007, col. 9: Ha narrow component flux, col. 10: Ha broad component flux. col. 11: Ha broad component 
FWHM, col. 12: [Nn]A6583 flux, col. 13: [Nn]A6583 FWHM, col. 14: optical Fell strength relative to HI3 broad component. 
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^ In units of erg s~ . 

* The nuclear light fraction is unreliable for FC < 0.25. 

Observed frame line flux in units of lO""*^^ erg cm 
In units of km . 
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Fig. 1.— SEE ATTACHED JPG FIGURE Representative examples of the EL-ICA starlight- 
nucleus decomposition. In each panel, we plot the observed spectrum (black), the combi- 
nation of the models (red), the decomposed components of the host galaxy (blue) and the 
power law continuum of the nucleus (green), and the Fell multiplets (pink). The observed 
spectra are smoothed with a boxcar of 5 pixels for the illustration. 

Fig. 2.— SEE ATTACHED JPG FIGURE Comparison of the narrow emission-line widths 
for ~ 3 000 high S/N type 2 AGNs. In the upper two panels. Ha, [Nil] and [SII] lines are 
all fitted with a single Gaussian. The profiles and redshifts of the [Nil] doublet are tied, 
and so are the [SII] doublet. The flux ratio of the [Nil] doublet is fixed to its theoretical 
value. It can be seen that FWHM(i7a) agree well with FWHM[NII] and FWHM[SII]. In 
the lower left panel, FWHM[NII] is compared to FWHM(H/3). The profiles and redshifts of 
the [Nil] doublet are tied to that of Ha and [SII] during the measurement of these lines. It 
can be seen that FWHM(H/3) agrees with FWHM[NII] within their respective uncertainties. 
FWHM(H/5) is plotted against FWHM[OIII] in the lower right panel. The scatter is much 
larger than those in the other three relations. 

Fig. 3.— SEE ATTACHED JPG FIGURE Sample results of the de-blending procedure 
applied to the H/3-|-[0III] (left panels) and Hq;-|-[NII] (right panels) ranges. 
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Fig. 4.— SEE ATTACHED JPG FIGURE Distribution of the relative measurement un- 
certainty of the H/9 and Ha broad component FWHM (upper left panel histograms). The 
correlation between the FWHM(Ha) and FWHM(H/5) of 1 220 NLSls with z < 0.4 in 
our sample are plotted in the upper right panel. The dash line denotes the best linear fit, 
FWHM{Ha) = 0.861FWHM{Hf3). The distribution of the discrepancy between the mea- 
sured and the expected values is plotted as the inset. The scatter of cj ~ 10% is only a little 
larger than the measurement uncertainty of FWHM(Ha) and FWHM(H/3) (5.5% and 2.6%, 
respectively). Our measurements of FWHM(H/5) are compared to those of WPM02 in the 
lower left panel. In the lower right panel, we plot the FWHM(H/3) ratio of our measurement 
to that of WPM02 against the height ratio of H/9 narrow to broad component. Objects 
with detected H/? narrow components are denoted by filled circles. Unfilled circles denote 
those in which narrow component of H/? is not detected. Note the concordance between our 
measurements and those in WPM02 for those objects in which the narrow component of B.(3 
emission-line is negligible. The difference between our measurements and those in WPM02 
in a handful of objects is due to the contribution of the H/3 narrow component, which was 
not considered by the WPM02 authors. 



Fig. 5.— SEE ATTACHED JPG FIGURE The upper panel shows the spectrum of SDSS 
J104755.93+073951.2, the only object in our NLSl sample with [OIII]/Hp > 3 and broad 
component of Ha and H/5 less than 2 000 km s~^. The broad permitted emission lines are 
illustrated on an expanded scale in the middle panel with the error array plotted in thin 
curve. The lower panel shows the newly identified type 2 SNe that fulfils our criterion of 
NLSls. It SDSS image is displayed in the inset panel. 



Fig. 6.— SEE ATTACHED JPG FIGURE Correlation between fiuxes of the Ha and H/3 
broad components. The dash line denotes the best linear fit to the data. Distribution of the 
Balmer decrement, the flux ratio of the Ha to H/3 broad component, is shown in the inset. 
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Fig. 7.— SEE ATTACHED JPG FIGURE Same as Figure 2 but for the 5 objects which 
were classified as NLSls in WPM02 but rejected as such in this paper. 

Fig. 8.— SEE ATTACHED JPG FIGURE Same as Figure 2 but for the 14 of the 86 objects 
classified in this paper as NLSl candidates but not included in the sample of WPM02. 

Fig. 9.— SEE ATTACHED JPG FIGURE Left panels: redshift and luminosity distribution 
of the NLSl sample (hatched area) and the parent BLAGN sample (open area). Right 
panels: dependence of the NLSl fraction on redshift and luminosity. 

Fig. 10.— SEE ATTACHED JPG FIGURE Luminosity versus FWHM for the broad com- 
ponents of the H/? and Ha emission lines, respectively. Note that the vast majority of NLSls 
populate in the region below the guiding lines (dashed) which represent a constant luminosity 
a few units of the Eddington luminosity. See §4.2 for detailed discussion. 

Fig. 11. — SEE ATTACHED JPG FIGURE Correlation between the nuclear monochromatic 
luminosity at 5100 Awith the H/? (left panel) luminosity and [OIII] (right panel) luminosity. 
The solid line is the best linear fit. 

Fig. 12.— SEE ATTACHED JPG FIGURE Distribution of the radio-loudness (left panels) 
and radio power (right panels) of NLSl and BLSls. 

Fig. 13.— SEE ATTACHED JPG FIGURE Dependence of the fraction of NLSls on the 
X-ray fiux (left panels) and X-ray luminosity (right panels). 

Fig. 14.— SEE ATTACHED JPG FIGURE Distribution of the relative strength of the Fell 
multiplets, -R4570, for the NLSl sample. Objects with 3 a upper limits are shown by the 
hatched histogram. 
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Fig. 15.— SEE ATTACHED JPG FIGURE Plots of i?457o against FWHM(H/3) (left panel) 
and FWHM(Hq;) (right panel). Objects with detected Fell emission are denoted as filled 
circles, and non-detections (3 a upper limits) are denoted as crosses. 



Fig. 16.— SEE ATTACHED JPG FIGURE Left panels: plots of H/3 luminosity versus the 
equivalent width of H/3 (upper) and FeII4570 (lower). Right panels show plots of monochro- 
matic nuclear luminosity versus the equivalent width of H/3 (upper) and FeII4570 (lower). 



Fig. 17.— SEE ATTACHED JPG FIGURE Plots of the soft X-ray photon index against 
the FWHM of H/3 (left panel) and Ha (right panel). 
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Fig. 18.— SEE ATTACHED JPG FIGURE The upper panel shows the distribution of 
the [OIII]X5007/HP^'^ ratio. Objects with rehably determined values are shown in the 
hallow histogram and those with undetectable flux of H/3 narrow component in hatched 
([O///]A5007/if/3"'^ = 10 is assumed for illustration purpose). Diagnostic diagram for the 
NLR of the 168 NLSls with reliably measurements of all the 4 narrow emission Unes are 
plotted in the lower panel. The solid line marks the theoretical classification line between 
Seyfert 2s and starforming galaxies (Kewley et al. 2000). 
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Fig. 19.— SEE ATTACHED JPG FIGURE Distribution of the relative redshift for [Oil] 
(a), [OIII] (b) and the broad components of H/3 (c) and Ha (d). The shift of [OH] is relative 
to the [Nil] and [SII] doublets. The shift of [OIII] and the broad components of Balmer 
emission hues is relative to the low- ionization forbidden hues, [Oil], [Nil], or [SII]. A few- 
tens of "blue outliers" are pinpointed in the inset. 

Fig. 20.— SEE ATTACHED JPG FIGURE Sample spectra of "blue outliers" (upper panel) 
and narrow line emitters (lower panel). Vertical lines indicate the expected centroid of 
emission lines according to the low-ionization forbidden lines. 

Fig. 21.— SEE ATTACHED JPG FIGURE Examples of "UV-deficient" NLSls. While most 
of these objects can be explained by intrinsic extinction indicated by large Balmer decrement 
(lower panel), the rest cannot (upper panel). 

Fig. 22.— SEE ATTACHED JPG FIGURE An example of "Fell-deficient" NLSls. Note the 
very narrow Balmer emission lines and also the very weak Fell multiplets. The prominent 
Hell broad emission hne ensures our classifying it as a type 1 AGN. 

Fig. 23.— SEE ATTACHED JPG FIGURE Spectra of 6 NLSls that possibly show 
the Mgll broad absorption trough. Three candidates, SDSS J102021.21+121909.1, SDSS 
J115816.72+132624.2, and SDSS J145724.01+452157.8, are selected from the SDSS DR4. 

Fig. 24.— SEE ATTACHED JPG FIGURE The relation between the stellar velocity dis- 
persion and the black hole mass for 308 NLSls for which a^, is reliably measured. The dash 
line is the relation given by Tremaine et al. (2002) for normal galaxies. 
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Fig. 25.— SEE ATTACHED JPG FIGURE Montages of 33 NLSls with rcdsliift less than 
0.1, for which either the host galaxies are viewed face-on or the SDSS 3" fiber aperture is 
dominated by the bulge contribution. 
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Fig. 26.— SEE ATTACHED JPG FIGURE Same as Figure 24 but for 33 nearby NLSls. 
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Fig. 27.— SEE ATTACHED JPG FIGURE The correlation between FWHM[NII] and a, 
for 206 NLSls with both parameters rehably measured. The sohd hne denotes the often 
used relation of FWHM[NII] = 2.35 a^. 

Fig. 28.— SEE ATTACHED JPG FIGURE Same as Figure 27 but for ~ 3 000 Scyfert 2 
galaxies with high S/N selected from the SDSS. The dash line marks the linear fit to the 
data. 

Fig. 29.— SEE ATTACHED JPG FIGURE Mbh-c^* relation for 613 NLSls in the sample. 
The dash hne marks the expected relation of Tremaine et al. using FWHM[NII] = 2.35 cr*. 

Fig. 30.— SEE ATTACHED JPG FIGURE Spectra of 4 NLSls with large amplitude vari- 
ability. Two epoch SDSS spectra are shown in red and green, and the difference in blue. 
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